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FOREWORD 

This document presents the results of work performed 
by personnel of the Fluid Mechanics Section of the Lockheed- 
Huntsville Research & Engineering Center, Huntsville, Alabama. 
This work was performed for the Aerodynamic Systems Analysis 
Section of the NASA-Johnson Space Center, Houston, Texas, (Con- 
tract NAS9-13429) in support of Space Shuttle plume impingement 
studies. The technical monitor for this contract is Mr. Barney 
B. Roberts, EX-32. 
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SUMMARY 

The Streamline Divergence Program was developed to demonstrate 
the capability to trace inviscid surface streamlines and to calculate "outflow 
corrected" laminar and turbulent convective heating rates on surfaces sub- 
jected to exhaust plume impingement. The analytical techniques used in 
formulating this program are discussed. A brief description of the Stream- 
line Divergence program is given along with a user's guide. The program 
input and output for a sample case are also presented. 
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Section 1 
INTRODUCTION 

Most current theories used in predicting convective heat transfer 
rates on body surfaces subjected to exhaust plume impingement do not 
utilize an adequate outflow correction theory. When a plume impinges on 
a surface at some angle of inclination, there is an outflow -induced thinning 
of the boundary layer, and a corresponding increase in the convective heat- 
ing rates to the surface. Application of current methods for predicting con- 
vective heat transfer (without an adequate outflow correction theory) often 
results in the predicted heating rates being considerably lower than the ex- 
perimentally measured heating rates. Without an outflow analysis to deter- 
mine the severity of the thinning of the boundary layer and the corresponding 
increase in convective heating rates it is very difficult to make experimental- 
theoretical comparisons. Under such circumstances, experimental -theoretical 
comparisons can result in erroneous conclusions. For example, convective 
heating rates are determined by both turbulent and laminar theory. The con- 
vective heating rates determined by turbulent theory are found to compare 
more favorably with experimental data than by laminar theory. The con- 
clusion is then made that the experimental data are turbulent which may' 
actually be laminar. When proper outflow correction theory is applied, such 
erroneous conclusions can hopefully be avoided. The results determined by 
laminar theory and by turbulent theory should compare favorably with the 
laminar and turbulent data, respectively. 

The computer program described in this report has been developed 
to: (1) demonstrate the capability to trace inviscid surface streamlines and 
calculate "outflow corrected" laminar and turbulent heating rates on a body 
subjected to rocket exhaust plume impingement (It is capable of analyzing 
impingement due to any supersonic flow field)j and (2) improve the existing 
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convective heat transfer prediction techniques (which do not apply outflow 
correction theory) of the Lockheed- Huntsville Plume Impingement Computer 
Program (PUMP) (Ref. 1-1). 

The analytical methods used in performing the inviscid surface stream- 
line and convective heat transfer calculations are described in Section 5 of 
Ref. 1-2. The Streamline Divergence program can be used to make: (1) com- 
parative evaluations of the effects of plume impingement due to different 
propulsion system configurations; (2) parametric studies; and (3) vehicle 
design studies. 
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Section 2 

PLUME IMPINGEMENT ANALYSIS USING THE STREAMLINE 

DIVERGENCE METHOD 

An analysis to trace inviscid surface streamlines and calculate '’out- 
flow corrected" heating rates on surfaces subjected to exhaust plume impinge- 
ment is a rather detailed procedure. The Streamline Divergence program 
is the means by which the final results are obtained. This section is intended 
to familiarize the program user with the data used by the Streamline Diver- 
gence program, and the means by which the data are generated. A brief com- 
putational scheme is presented in the following paragraphs. 

The prominent features of the computational scheme required in order 
to completely describe the environment of a vehicle subjected to exhaust plume 
impingement, and the effects of this environment on the vehicle are: 

• Description of the thermochemical behavior of the 
propellant system, 

• Prediction of the undisturbed exhaust plume properties 
via gas dynamic and thermochemical methods , 

• Prediction of the body local flow properties for a 
surface immersed in an exhaust plume, and 

• Prediction of heating rates to the surface. 

The thermochemical properties are calculated using the NASA-Lewis 
TRAN72 computer program (Ref. 2-1). The free, undisturbed plume local 
flow properties are computed using the method of characteristics (Ref, 2-2), 
with the computations accomplished using the Lockheed-Huntsville Variable 
O/F Method-of-Characteristics (VOFMOC) computer program (Ref. 2-3). 

The body local flow properties for a surface immersed in a free plume are 
calculated using the PLIMP program. Finally the heating rates to the sur- 
face are calculated using the surface geometry and impingement properties. 
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2.1 THERMOCHEMICAL, ANALYSIS 

The thermochemical properties of the propellant system being analyzed 
are computed assuming the combustion products to be chemically frozen or 
in local chemical equilibrium. The thermochemical properties are obtained 
by computing the flowfield expansion properties from specified engine com- 
bustion chamber conditions using the NASA-Lewis TRAN72 computer program. 
The TRAN72 program has been modified to generate a binary tape containing 
all the pertinent parameters needed for flowfield and heat transfer analysis 
including the mole fractions of each constituent. This method reduces much 
of the tedious and time-consuming effort normally necessary to prepare input 
for a single flowfield calculation. 

2.2 UNDISTURBED FLOWFIELD ANALYSIS 

The VOFMOC program is used to calculate the exhaust plume flow fields. 
Utilizing the VOFMOC program, the exhaust plume for any arbitrary nozzle 
and exhaust gas (ideal or reacting) system with gradients in the oxidizer/fuel 
ratio can be calculated. At finite back pressures the program normally gene- 
rates the plume flow field using the flow regime continuum methods. As the 
back pressure is lowered, the outer regions of the plume become progres- 
sively more rarefield until eventually the flow becomes free molecular. The 
VOFMOC program permits transfer of the calculations from a continuum 
analysis to a free molecular analysis across a "freeze" front in the flow field. 

2.3 BODY LOCAL FLOW PROPERTIES ANALYSIS 

The PLIMP program is used to calculate body local flow properties for 
a surface immersed in a free plume. The PLIMP program locates each im- 
pinged body and the particular motor with respect to a reference coordinate 
system. The impinged bodies are then further divided into small elemental 
areas over which uniform approach flow is assumed. The geometrical sub- 
shapes which the computer program can analyze are shown on the following 
page. 
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• Conic shapes, i. e., cones, cone frustums, cylinders 

• Flat plates, i. e., rectangular, circular, triangular 

• Airfoil shapes, and 

• Arbitrary axisymmetric bodies whose surface curva- 
ture can be described by a polynomial curve fit 
(maximum of fourth degree). 

The local flow conditions on a given elemental area are obtained by 
locating the centroid of each elemental area in the free plume from a search 
of the plume fLowfield tape. The local impingement angles are obtained using 
the local flow velocity vector and the elemental area unit-normal vector. 
Shock calculations are then performed to obtain the local flow properties. 

An iterative solution is used to perform the equilibrium shock calculations 
for a real or ideal gas. The real and ideal gas calculations are similar, the 
difference being that the ideal gas solution converges on the first iteration. 

A more detailed description of the equilibrium shock calculations is given 
in Ref. 2.4. Forces and moments are then obtained by numerically integrating 
the local impact pressure distribution over the body. 

Pressures can be calculated in continuum flow using Newtonian theory; 
modified Newtonian theory; or oblique shock theory; in free molecular flow 
using kinetic theory; and in transitional flow using empirical relationships. 

Convective heating rates in the continuum or transitional flow regimes 
are calculated (outflow correction theory is not applied). Heating values for 
continuum-laminar, continuum-turbulent and transitional flow are calculated 
The method used for the calculation of heat transfer to a body in the trans- 
itional regions of the plume has been previously presented in Ref. 2.4. Basic- 
ally, it amounts to a continuum solution correction to account for the effects 
of a slip velocity and temperature jump at the surface. Free molecular heat- 
ing rates are calculated at each elemental area from an energy balance- type 
equation. 
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ZA STREAMLINE DIVERGENCE ANALYSIS PROGRAMS 

The computer programs and calculational techniques described pre- 
viously have been combined in an automated fashion to fully describe the 
vehicle environment and the effects of this environment on the vehicle. A 
schematic flow chart of the calculational process and the flow of information 
is given in Fig. 2- 1. The arrangement in Fig. 2- I allows the greatest flexi- 
bility with a minimum of effort and prevents redundant calculations. For 
instance, if the heating effects are desired for a previously developed plume, 
but a new body, then only the PLIMP and Streamline Divergence programs 
need be rerun. 
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Fig. 2-1 - Calculational Process 


2-5 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 





LMSC-HREC TM D390450 


Section 3 

CALCULATION OF CONVECTIVE HEATING RATES USING 
A STREAMLINE DIVERGENCE METHOD 

This section describes the analytical procedures used by the Stream- 
line Divergence program. The analysis permits the tracing of inviscid 
surface streamlines, and the calculation of "outflow corrected" convective 
heating rates on surfaces subjected to exhaust plume impingement. Basic- 
ally, the procedure is the following: the Streamline Divergence program 

applies the Method of the Bicubic Piecewise Polynomial Functions to the 
impingement pressure distributions calculated by the PLIMP program so 
that the pressure distribution over the entire body may be described by 
functional means. The same calculations are applied to the local velocity, 
local entropy, freestream velocity and freestream entropy distributions 
determined by the PLIMP program. Once these calculations are performed, 
the body local flow properties are defined over the entire body in functional 
terms. The Streamline Divergence program next applies DeJarnette' s {Ref. 
3-1) streamline divergence methods to determine the streamline patterns 
and corresponding "equivalent radius" or metric coefficient for the particular 
body surface subjected to exhaust plume impingement. Convective heating 
rates for a laminar and turbulent boundary layer are calculated using the 
integral form of the energy equation. The effects of variable freestream 
velocity, density and pressure are accounted for through the use of the ap- 
propriate transforms of the flat plate solution. Reynolds number is output 
in order for the user to choose the most applicable of the two convective 
heating rates that are provided. 

3.1 METHOD OF THE BICUBIC PIECEWISE POLYNOMIAL FUNCTIONS 

The tracing of inviscid surface streamlines and calculation of convective 
heating rates are dependent on the surface pressure distribution. Before in- 
viscid surface streamlines and heating rates can be calculated, a description 
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of the pressure distribution on the body surface is needed. This pressure 
distribution must yield continuous first and second derivatives in both the 
axial and circumferential directions (assuming that a cylindrical coordinate 
system is used). If the pressure distribution on the body surface could be 
described by an analytical equation this requirement could easily be met. 
Unfortunately, this is not the case; most if not all pressure distributions 
cannot be described by an analytical expression. Often, only the pressure 
at several points along the circumference of the body surface at several axial 
stations is known, and the pressure distribution must be generated from these 
data. Therefore, a numerical technique is required that will generate the 
pressure distribution from these data, and calculate the required derivatives 
at "any" point on the body. This technique must also be capable of performing 
numerical interpolation and differentiation with a reasonable degree of accuracy. 

Experience has shown that numerical differentiation methods are gener- 
ally inaccurate. The method of splines (Ref. 3-2) however, has proven to be 
an effective method for numerical interpolation, differentiation and integration. 
The two-dimensional or bicubic spline (which is needed for a three-dimensional 
axi symmetric body) is able to perform these operations accurately due to the 
strong convergence properties it possesses. 


The Streamline Divergence program applies the Method of the Bicubic 
Piecewise Polynomial Functions to the pressure distribution determined by 
the PLrIMP program so that the pressure distribution may be described in a 
functional means. The governing equations, as derived in Ref, 3-2 may be 
written as follows: 

= c oo(;s) (a?) +c io(s) (s?) +c oi(s;) (^) 

L=0 j=0 
3-2 


(3.1) 
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E E (£)' 

i=2 j =0 


(3.6) 


where: 


P = local flowfield property of interest 
C i - = coefficients for the polynonial P (x, <jt) 

J 


and the corresponding coordinate system used is presented in Fig. 3-1. The 
same equations are applied to the local velocity, local entropy, freestream 
velocity and freestream entropy distributions determined by the PLIMP pro- 
gram. Once these calculations are performed, the body local flow properties 
are defined over the entire body in functional terms. The end result obtained 
from applying the Method of the Bicubic Piecewise Polynomial Functions to 
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Fig. 3-1 - Coordinate System for Pressure vs X and <f> 
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the data presented in Fig. 3-1 is presented in Fig. 3-2. The pressure surface 
is completely specified by individual sets of coefficients which uniquely 
describe a bicubic polynomial function in each region over the entire pressure 
surface. At any point on the surface, the pressure and pressure derivatives 
can be determined by using Eqs. (3.1) through (3.6) and the corresponding set 
of coefficients. 


Briefly, to determine the value of a particular flowfield property at a 
point defined by X and <j> on the body surface, the region which contains 
x and <j> is first specified. Equations (3.1) through (3.4) are then evaluated 
at each of the four corners of the region R„. These equations make up a 
set of sixteen simultaneous equations to be solved for the sixteen coefficients 
describing the particular flowfield property of interest in region R... The 16 
coefficients are solved for, and then with the value of X and <t> substituted into 
Eq. (3.1) to yield the particular flowfield property value. 


3.2 THE STREAMLINE DIVERGENCE METHOD 


At this point, the body local flowfield properties are defined over the 
entire body in functional terms. The Streamline Divergence program next 
applies the methods of Dejarnette (Ref. 3-1) to determine the streamline 
patterns and calculate the corresponding "equivalent radii" or metric coef- 
ficients along each streamline for the particular geometry subjected to ex- 
haust plume impingement. The general form of the equations, as derived 
in Ref. 3-1 may be written as follows: 


D§ 

DS 


pv‘ 


sin0 cosT ^ — 
ox 




(cosG cos6 . + sin0 sin6 , sinT) 

+ £ _ £ 1 

f d<J> 


sinF 


COS0 cosr -5“ + 
ox 


(sin0 cos6^ - cosG sin5 , sinr) 

7 2 af 


(3.7) 
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Fig. 3-2 - Smooth Pressure Surface Fit Generated by the Method of the Bicubic 
Piecewise Polynomial Functions 
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where 


DS 

p 

s 

P 

V 

9 

P 

f 

x 

<f> 

g 

r 


differential arc length along a streamline 
stagnation pressure 


mass density 

inviscid speed on surface 

inclination angle of the inviscid surface streamline 
static pressure 

body radius, measured from the longitudinal axis 
longitudinal coordinate 
circumferential angle 



body angle defined by sinF = 


3_f 
B x 
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1 8f 

T' 

° - <t> ~ 6 <t> 

0 = coordinate normal to the streamline on the body surface 

h = "equivalent radius" or metric coefficient in the |3 direction 
M = Mach number 



6^ - body angle defined by sin6 ^ = 




1 + -r 



For a cylinder, which is the only geometry currently considered by the Stream- 
line Divergence program; the angles T and 6^ are identically equal to zero, 
reducing Eqs. (3.7) through (3.10) to: 
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respectively. It should be noted at this point that other geometries can 
easily be incorporated into the Streamline Divergence program simply by 
making the proper substitutions into Eqs. (3.7) through (3.10) to yield new 
applicable Eqs. (3.11) through (3.14). 

The values of local static pressure, velocity and the derivatives 

_9p ap j^p a z p d 2 p 

d * ’ 3x 2 ' 80 2 w 

are determined by using the bicubic piecewise polynomial functions at the 
correct geometric location on the body surface. The inviscid surface stream* 
lines and their corresponding metric coefficient, h, are calculated by numer- 
ically integrating Eqs. (3.11) through (3.14) simultaneously. The numerical 
integration scheme used is the fourth order variable step size Runge-Kutta 
method. The initial conditions required to start the integration of each 
streamline are determined external to the program and input by means of 
punched data cards. 

3.3 RELATIONSHIP OF THE METRIC COEFFICIENTS TO CONVECTIVE 
HEAT TRANSFER THEORY 

The convective heating rates for a laminar and turbulent boundary layer 
are calculated using the integral form of the energy equation. The effects of 
variable freestream velocity, density and pressure are accounted for through 
the use of the appropriate transforms of the flat plate solution. Non-constant 
properties through the boundary layer are also accounted for. 

In Ref. 3-3 (among others) it is shown that the energy equation can be 
reduced to an equation with enthalpy as the dependent variable under a number 
of different conditions including some cases where there is dissociation, 
chemical reaction and mass diffusion within the boundary layer. The com- 
plete development of the enthalpy form of the equation and the evaluation of 
the surface heat rate are presented in Ref. 3-3, 
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Based on examination of several ''exact' 1 laminar boundary layer solu- 
tions, Eckert (Ref. 3-4) recommends that the effects of variable gas properties 
through the boundary layer can be accounted for by simply evaluating the prop- 
erties at a "reference enthalpy" and using these values in the constant property 
solutions as obtained by Blasius. Based on this method, the convective heating 
rate to the wall is evaluated using: 


Laminar Boundary Layer 



Pr 


0.332 
* 2/3 


(P* M* V e ) 

^ -5 


0.5 


{H -h ) 
1 r w 


(3.15) 


where: 


laminar characteristic running length. is obtained for 

variable property flow over an arbitrary shape, by numerically 
integrating the following equation along a flowfield streamline, 
(Ref. 3-3) 



* 

P 


* 2 
P V h 



V h 2 ds 
e 


= adiabatic wall enthalpy 


H 


r 




r) 


r = recovery factor 


r 
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Turbulent Boundary Layer 


q = 


where: X^, 


Re 


T 


r 



. . (P* (il)- 2 

1 + .212 fReJ • (P r - 1) 6 ' X T / 


(H -h ) 
r w' 


= turbulent characteristic running length 


P u t v e h 


hrs f p%*v e h l - 25 d S 


turbulent Reynolds number 

Re T = p V e X T /jU* 

p * 1/3 

r 

Prandtl number 
absolute viscosity 

velocity at the edge of the boundary layer 

freestream velocity 

wall enthalpy 

freestream enthalpy 

wetted length along the streamline 

metric coefficient 


(3.16) 


The starred (*) properties refer to properties evaluated at a temperature 
corresponding to a "reference enthalypy," H*, where 
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H* = h + 0. 5 (H - h ) + 0.22 (H - h ) 
e o e' 're' 

h g = static enthalpy at the boundary layer edge 
H q = stagnation enthalpy 

The wall enthalpy (h w ) in Eqs. (3.15) and (3.16) may be calculated 
assuming either an infinitely fast recombination rate in the boundary layer 
(chemical equilibrium) or an infinitely slow recombination rate (chemically 
frozen). The upper extreme to the convective heating rate is determined 
by assuming equilibrium chemistry throughout the boundary layer and a 
fully catalytic wall. The lower extreme to the convective heating rate is 
determined by assuming a frozen thermal boundary layer and a non- catalytic 
wall. 


Eckert’s reference enthalpy method has been chosen for use in the 
above analysis for demonstration purposes. Other methods do exist (though 
they have not been evaluated under the present study) that may be applied to 
the problem of exhaust plume impingement and, can be incorporated into the 
computer code with a minimum of programming difficulty. 
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Section 4 

USER'S INPUT GUIDE FOR THE STREAMLINE 
DIVERGENCE COMPUTER PROGRAM 

The Streamline Divergence program utilizes three forms of input: (1) 
data statements within the program; (2) magnetic tapes; and (3) punched data 
cards. Data statements are used to input thermodynamic data (for reaction 
products in functional form) that seldom change. Magnetic tapes and punched 
data cards are used to input data that frequently vary from one exhaust plume 
flow field to another or from run to run. The program magnetic tape assign- 
ments are given in Table 4-1. Each punched data card and its use is explained 
in Section 4. 1 . 

4.1 PROGRAM INPUT GUIDE 

Card 1: Program Control Card 


Item 

NFIT, number of flowfield properties to be 
curve fit using the Method of Bicubic Piece- 
wise Polynomial Functions (upper limit of 5) 

IUNIT, tape unit number on which impinge- 
ment data are written 

NOXCUT, number of axial printout intervals 

NOYCUT, number of angular printout intervals 

JMAX, number of points on each streamline 
where heating rates are to be calculated 
(upper limit of 150) 

ICARD, number of species data input on 
punched data cards 

NPRINT, control option that indicates to print 
(= 1) or not to print (=0) intermediate data 
(used during program checkout procedures) 


Format Column 

1615 5 

10 

15 

20 

25 

30 

35 
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Table 4- 1 

MAGNETIC TAPE ASSIGNMENTS FOR THE STREAMLINE 
DIVERGENCE PROGRAM 


Where Required 

Tape Units 
U-1108 

Tape Unit Function 

Section 1 

Subroutine PLIRED 

Subroutine PLIRED is 
the controlling routine which 
arranges the PLIMP output 
in the form used by the 
data acquisition routines. 
This subroutine is called 
once at the beginning of 
each run. 

IUNIT 

Impingement data generated 
by the PLIMP program (input 
data) 

Section Z 

Convective Heating Rate 
Analysis 

This section performs the 
convective heating rate 
analysis along each inviscid 
surface streamline being 
traced. 

3 

Thermochemical properties 
data generated by the TRAN7Z 
program (input data) 
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Card 2: 

F ormat 
1615 

Card 3: 
6F10.4 


Card 4: 
10E11.4 

Card 5; 
1615 

Card(s) 6 

4{A6, 4X, 

El 0.6) 


Column 


Item 


5 NOSE, number of individual streamlines to 

be traced {upper limit of 10) 

This card contains the initial conditions required to start the 
tracing of each streamline utilizing the fourth order variable 
step- size Runge-Kutta method. There is one card for each 
streamline to be traced. 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 


XSTART, initial axial location of streamline (in. ) 

YSTART, initial angular location of streamline 
(rad) 


TSTART, initial streamline direction (rad) 

HSTART, initial metric coefficient (in.) 

ASTART, initial derivative of metric coefficient 
with respect to streamline coordinate (in./in.) 

DELS, integration step size along streamline (in.) 


TWALL, body surface wall temperature. Assumed 
to be constant over the entire body (°R) 


5 NENTLP, control option that indicates to calculate 

wall enthalpy assuming frozen (=1) or equilibrium 
(=2) chemistry throughout the boundary layer. 

■: Species concentration cards (used when NENTLP = 2, and 
1CARD > 0). This card(s) specifies the mole fractions of the 
flowfield gas composition. 

1"6 Chemical species name, i.e., H-,, H„0 

21-26 2 2 
41-46 
61-66 

11 “20 Mole fraction 

31-40 

51-60 

71-80 
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The species names are left-adjusted. 


Cards 7, 8 and 9: Wall enthalpy array cards (used when NENTLP = 1). These 

cards are used to define wall enthalpy as a function of impinge- 
ment pressure and body surface wall temperature. 


Card 7: 


Format Column 
515 5 

10 


Card(s) 8 and 9: 


Item 


NPWALL, number of independent pressures in 
the impingement pressure array 

NTWALL, number of independent temperatures 
in the body surface wall temperature array 


E10.6/ 

(8E10.6) 


Card 8: 


TWALL, temperature in body surface wall 
temperature array (°R) 

Card(s) 9: 


PWALL, pressure in the impingement pressure 
21-30 array (psi) 

41-50 

61-70 

11-20 HWALL, corresponding dependent wall enthalpy 

31-40 (Btu/lbm) 

51-60 

71-80 


Cards 8 and 9 are repeated NTWALL times. 
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Card 10: This card contains identifying information for the particular problem 

of interest. This information is printed at the top of each page of 
the output. 


Format 

Column 

Item 

12A6,I5 

1-72 

Comment card; header information such as 
problem title may go in these columns. 


73-77 

Case number of the particular problem being 
analyzed. 
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Section 5 

PROGRAM NOMENCLATURE 

The Streamline Divergence Program, like all computer programs, 
contains its own particular system of terms. This section presents a list 
of all "key" nomenclature used in the program, and may be used as a quick 
reference when working with the program. 


5-1 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TM D390450 


A : , i a i i 


DAD-> 


, u E 1_ ; j 


D H L j 


■u ruo 


.DXus 
UYlKj 
A A if- 
LAW! 

L W 

HS T AR 1 

HE 

HWALL 

H2 

I CArfU 
1 UN I T 
JMAX 

NOlA 

NEN T LP 

NF I I 

NOSL 
NO XU UT 
NOYUUT 


INITIAL FI Kb I uiPIVAlIVE OF ME I R 1 C COEFFICIENT WITH RESPECT TO 
OISTAnCl ALONG A ST PL AMU 1 Nt (REQUIRED TO STAPT THE TP AC I NG OF EACH 

streamline ) 

bLCOND DERIVATIVE OF ME TP I C COEFFICIENT WITH PEbPECT TO DISTANCE 
ALONG A STREAMLINE 

INTEGRA T ION STEP SIZE ALONG A STREAMLINE (INCHES) 

FIRST DcRiVATIVi.:. OF METRIC COEFFICIENT WITH RESPECT TO DIbTANCE 
ALONG A STRLAMLINE 

FIRST DERIVATIVE OF STRLAMLINE DIRECTION WITH RESPECT TO DISTANCE 
AL ONG AST Rl AML 1 NL 

FIRST l)lKWaT[Vc. OF X w IT H RESPECT To DISTANCE ALONG A STREAMLINE 
FIRST ULR1VATIVE OF PHI WITH RESPECT TO OiSTANCE ALONG A STREAMLINE 
LAMINAR ADIABATIC WALL ENTHALPY (BTU/LtJM) 

TUPHULliWT ADIABATIC WALL ENTHALPY (BTU/LbM) 

WALL ENTHALPY (BTU/LBM) 

INITIAL Me I PIC COEFFICIENT OF STREAMLINE (REOUIRED To START THE 
TRACING OF EACH STREAMLINE) (INCHES) 

ENTHALPY AT THE LOGo OF THE BOUNDARY LAYER (BTU/LBM) 

ENTHALPY I ri THE WALL ENTHALPY ARRAY (bTu/LBM) 

METRIC uOlF-F I C I ENT (INCHES) 

NUMBER OF SPECIE DATA INPUT ON PUNCHED DATA CARDS 

tape UNIT NUMBER ON WHICH IMPINGEMENT DATA IS WRITTEN 

NUMB UR OF POINTS ON EACH STREAMLINE WHERE HEATING RAT lS ARE TO BE 
CAECULA f ED 

COUNTER That INDICATES the NUMBER OF The INVISCID SURFACE STREAMLINE 
PRESENTLY BEING TRACED 

control option that indicates to calculate, wall enthalpy assuming 
frozen (=1) op EQUILIBRIUM (=2) CHEMISTRY THRUuGHOuT THE BOUNDARY 
LAYER 

NUMBER OF FLOW FIELu PROPERTIES TO BE CURVE FIT USING THE METHOD 
OF BICUBIC PIECEWISE POLYNOMIAL FUNCTIONS 

NUMBER UF INDIVIDUAL STREAMLINES TO BE TRACED 

NUMBER Up AXIAL PRINT OUT INTERVALS 

NUMBER UF ANGULAR PRINT OUT INTERVALS 


original pmseis 

OF POOR QUAUTS 
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NPH J NT 
NP w^LL 

n r wall 

PHI 
PR 
PR S 
PRST 
PW 

P Vi al.L 
<jl 

OT 

PE LAM 
k L TUPLJ 
RHU 
R Hub 

I 

5 

S I Nh 

It 

TK 

T^l ART 
TW 

TWALL 

U 

V ( A U) « C ’ D 

V ( - 

V ( - ,- »[) 
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CONTROL OPTION THAT INDICATES TO PRINT (=1) OR NOT TO PRINT (=0) 
INTERMEDIATE data 

NUMBER uF INDEPENDENT PRESSURES IN THE IMPINGEMENT PRESSURE ARRAY 

NUMBER of independent temperatures IN The body surface wall 
TEMPERA TU KE ARRAY 

ANGULAR LOCATION ON THE BODY SURFACE (RADIANS) 

LOCAL PKAROTL NUMBER 

REFERENCE LAMINAR LOCAL PRANDTL NUMBER 

reference Turbulent local prandtl number 
lucal impingement pressure (psf> 

PRLSSuPL IN the impingement PRESSURE ARRAY (PS I) 

LAM I NAP CONVECTIVE HEATING RATE ( BTU/FT**2 SEC) 

TURBULENT CONVECTIVE HEATING RATE (bTU/TT**2 SEC) 

LAM I NAP REYNOLDS NUMBER 
TURBULENT REYNOLDS NUMBER 
LOCAL DENSITY ( LBM/FT**3) 

REFERENCE LAMINAR LUCAL DENSITY (L.BM/FT**3) 
reference Turbulent local ulnsity (lbm/ft**3) 

LOCAL ENTROPY (BTU/LBM*R) 

FREE STREAM ENTROPY (BTU/LBM °R) 

LOCAL FLOW FIELD TEMPERATURE AT THE EDGE OF THE BOUNDARY LAYER ( °R ) 
BUOY SURFACE WALL TEMPERATURE ( °K ) 

INITIAL STKeAMLINE DIRECT I UN (REUUlRLU TO START THE TRACING OF 
EACH STREAMLINE) (RADIANS) 

BODY SURFACE WALL TEMPERATURL ( °R ) 

BODY SURFACE WALL TEMPERATURE. ASSUMED TO BE CONSTANT OVER THE 
ENTIRE BODY (°R) 

LOCAL VELOCITY (FT/SEC) 

> ARRAY UotD IO STORE LOCAi_ Fi_uW PROPERTIES ANU THEIR DERIVATIVES AS fi 
FUNCTION UF GEOMETRIC LOCATION ON THL BODY SURFACE 

) DEFINES the AXIAL GeOMl TRIG LOCATION ON THE BODY SURFACE (1NCHLS) 

) DLF“lNL-> The ANGULAR GEOMETRIC LOCATION ON The BODY SURFACE (RADIANS) 

ORIGINAL PAGES 5-3 
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V(-tBt-)-) 


V(A« — « — »_) 


V I S 
VINF 

V I ss 
VlSST 
VMAX ( 1 ) 
VMaX ( 2 ) 

V 1*1 AX ( 3 ) 
VMAX ( 4 ) 
VMAX ( 5) 
X 

XLLAW 

XL TUPB 
XM 

XSTART 
XX 1 

YDUT ( 1 ) 
YDUT ( 2 ) 
YDUT C 3 ) 
YD OT ( 4 ) 
YDUT (5) 
YS1 AR I 

V Y 1 


DEFINES THE DERIVATIVE UF V WITH RESPECT TO X (B=2>» WITH RESPECT 
TO PHI (B=3)» OP WITH RESPECT TO X AND PHI <B=4>* IF B =]< NO 
DERIVATIVE HAS BEEN TAKEN 

defines Ihc local flow, ppopekty, local velocity <a=i>*lqcal entropy 

<A = 2). FREE ST PE AMI VELOCITY <A=3)» FREE STREAM ENTROPY (A = 4), AND 
COcAL IMPINGEMENT PRESSURE (A = S) 

AisSOluTL VISCOSITY (LbM/FT SEC) 

FPlE STREAM VELOCITY (FT/SLC) 

REFERENCE LAMIINAP ABSOLUTE VISCOSITY (EBM/FT SEC) 

reference Turbulent absolute viscosity < lbm/ft sec) 

MAXIMUM LuCAL VELOCITY UN THE BODY SURFACE (FT/SEC) 

MAXIMUM LOCAL ENTROPY ON THE BODY SURFACE < BTU/LBM °R) 

MAXIMUM FKLE STREAM VELOCITY (FT/SEC) 

MAXIMUM FKLE STREAM ENTROPY (BTU/LBM °R ) 

MAXIMUM LOCAL IMPINGEMENT PRESSURE ON THE BODY SURFACE (PSF) 

AXIAL LOCATION ON THL BODY SURFACE (INLHLS) 

LAMINAR CHARACTLRI ST 1C RUNNING LLNGTH (INCHES) 

TURBULENT CHARACTERISTIC RUNNING LENGTH (INCHES) 

LOCAL MACH NUMBER 

INITIAL AXIAL LOCATION UF" STREAMLINE (REUUIRED TO START THE TRACING 
OF EACH STREAMLINE) (INCHES) 

AXIAL LOCATION ON THE BODY SURFACE (INCHES) 

thl_ stoked value of dxds 

THc STOKED VALUE OF DYDS 

the stored value of dtds 
TH c. STORED value OF DADS 
THE STOKED VALUE OF DHDS 

INITIAL ANGULAR LOCATION of STREAMLINE (REUUIRED TO START THE 
TRACING OF EACH STREAMLINE) (RADIANS) 

ANGULAR LOCATION ON THE BODY SURFACE (RADIANS) 


ORIGINAL PAGE IS 
OF POOR QUALITY! 
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Section 6 

USE OF THE STREAMLINE DIVERGENCE PROGRAM 

To familiarize the user with the operation of the program, a sample 
case is presented. The user is guided through a step-by-step procedure 
in setting up the data, coding the input for the program and interpreting the 
output from the program. 

6.1 SAMPLE PROBLEM 

A typical shuttle configuration is chosen as the example. In this ex- 
ample, one is interested in tracing inviscid surface streamlines, and calcu- 
lating "outflow corrected" convective heating rates on that portion of the 
booster fuselage subjected to exhaust plume impingement. The overall coor- 
dinate system and the orientation between the booster and the orbiter engine 
is shown in Fig. 6-1. The coordinate system reference point corresponds to 
the first point analyzed by the PLIMP program. The upper surface of the 
booster fuselage was mathematically modeled as a cylinder (radius of 5.805 
in.). 

6.2 IDENTIFICATION OF THE INPUT/OUTPUT DATA 

The input data for the problem configuration described in Section 6. 1 
will be illustrated by identifying the computer printout of the input data in 
Table 6-1. 

The input data for the configuration of Section 6.1 is identified as 
follows: 

© This is an identification of the problem, and is printed at 

the top of each page of data. This information is contained 
in card 10 of the input. 
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Table 6-1 

STREAMLINE DIVERGENCE PROGRAM OUTPUT LISTING 


PAGE 1 



■ ;**' 


CASE 

•NO • 1 





(l) LOCKHEED /HUNTSV j LLF 

STREAMLINE 

D I VERGEnCE 

COMPUTER 

program Sample 

OUTpUT 

( 

2) NF IT 

1UN 1 T 

NOXCUT 

NOyCUT 

jmax 




5 

1 l 

20 

2S 

1 4 S 




r STREAMLINE 

X ST ART 

YSTART 

TSTART 

hstarT 

astart 

DELS 


I 

3.0603 

.ocoi 

.050 3 

5.8050 

.000c 

.1030 

2 

3 • C 6Q0 


• 1030 

5 . ROSO 

. 000 3 

.1000 

0 s 

3 

3.06C3 

*0005 

* » sot- 

5.8050 

.0000 

.1030 


4 

3.06C3 

* 0QC3 

.20 30 

5.B050 

.0030 

. 1 000 


• S 

3*0601 

.3000 

» 2S0Q 

5 . 8 © 5 it 

.0030 

. 1C30 


TWALL NENTlP 1CAHD 

(?) 540.0000 2 C 

TOTAL FREESTREaM ENTHALPY = -.8667+OC BTU/LBM 
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Table 6-1 (Continued) 



O' 

I 


ST«t*»iS*C 




A 

I INI 
Pi 

|f$f I 


CASE VO. | 

LOCKMCeo/rtUNTSVlLlC STlEA«LlNF OlVERlErjCE 

Phi velocity ENTROPY nacm 

|p*0l IfT/StCJ <FT**2/SEC*+2/P> 


computer prosmam sample output 
mar 


RHO 

ILAM/Cf l< 


V I S A PR A 

4 L 8 M/F T/SET 4 


ERUIv R 


XLTURB 


/^\ r^\ 

|7R«-02k§y* l232 + ?5v3'.2C39*''5 ' v i9i 29 0 | + « 


•3159+U .3S!|-?2 

•99*5+01 .19*6-01 


•1262+15 

• 56 * 1 +*^ 


* 58,T * c, ^ , * ,, ^ e (2 


• 187 6*15 

• S*53 + *M 





• 271* + IJ0 


•J4S9+3I 

•5374-12 

■!209*}S 

.172^+05 

• 322 S +*| 

• 2100-01 

+9016+^| 

•1851-31 

♦ 5729 + 0-1 

•S9J 1 + 0 | 

. ISRl+oO 

•3591+OC 

*>•59+31 

•7218-?? 

*1313+05 

• 1S91+05 

. 3339+o 1 

• 2507-01 

*tot«+o 2 

•1702-01 

•S7*6+?v 

•5901+0] 

•1215+-C 

•1291+00 


* 9 1 55-52 

•1315+05 

• 1175+25 

•3157+01 

•2919-31 

+ 1121+12 

.1523-01 

♦5815+00 

♦+902+01 

•15S1+-3 

•1929+00 

• 12 S 6 +?| 

•1109-ci 

• 1351 + 05 

•1372+35 

•3547+rl 

•3307-01 

•1213+12 

♦136S-0H 

•5077+00 

# 6 1 9o*h | 

•5572+09 

•5552+00 

•115*+n 

•1307-01 

•IJ’l+OS 

* 1 27 + + *5 

. 3 * 68 +— | 

•3442-01 

• 1 295+12 

•1228-01 

• 59* l+ ? r. 

. 6310*0 t 

•4115+^3 

• 6 | 02 + CO 

.9450+31 

•1507-01 

• 139S + 05 

• i 1 9fl*05 

• 3762+M 

•3901-01 

* l 364+12 

• 1108-01 

. 59+7*o« 

• 411 2 + h | 

♦6723+09 

• 6 * 2 1 +00 

*l#S7+il 

• 1 7 1 0*0 1 

•1396+05 

• 1 127 + 25 

• 3019 + M 

•1271-01 

•1927+12 

•1902-nl 

• 5975 + 00 

.6501+01 

. 7123 + + 1 

.7589+00 

• 50S7 + H 

•1914-31 

• 1 1^9*05 

• 1C67 + 05 

.3931+01 

.1530-01 

• 1*70 + 12 

.3907-01 

♦4071+00 

•6735+01 

.7973+^0 

.8213+00 

• S257+H 

• 2 | 2 l -01 

• 1 vJr+jS 

• I0C2+05 

•1316+ * | 

•1781-01 

•1*19+1* 

• 3*29-?l 

• 6321+00 

• 68 9, 1 + 01 

■*$p?+~p 

•0021+Or 

• 5154+i| 

•2335-01 

•1133+05 

.9175+01 

• 10*3+M 

•sroi-oi 

•1552+32 

• 371 1— ?! 

• 4316+10 

• 7 3 6 1 ♦ *!* 1 

•9325+-0 

•9125+0? 

•5*54+4) 

• 2519-c l 

•I130+7S 

.8947+01 

• 1 1 19 + - | 

• 52,75-1?! 

•1S77+-2 

•3460-? 1 

♦6366+00 

.7234+01 

• 9537* - 0 

• IP l2 + n l 

• 5*55++| 

• 2765— C 1 

■ 1 l*7*o5 

• 0199 + 41 

.1216 + M 

•5331-01 

♦1593+12 

• 3439- 1 :! 

•60*l+?0 

*7116+01 

• i?:5+m 

• 1 ?6 3+0 1 

♦ 4955+ H 

•29*1-01 

♦ IH51+R5 

•0068+ri 

.120-3+-I 

•5532-Cl 

♦ l*02+*2 

.3538-01 

• 61 M+Ot 

.7603+01 

• 1 ?5 7 ♦ - 1 

• 1 121+01 

■ «25S+ - | 

.3205-" 1 

« H6I+05 

.767 t +01 

• 1 31 2+ * | 

•5651-Pi 

• l*ll+*2 

.3!79-;s 

• 6 1 t 6+ 0 r 

,7797 +^ 1 

. 1 l^9+^| 

♦ 1 187 + ^1 


*• 7581*03 

. 2641 * 0 ? 

- . 78 1 6+q3 
• 2 * 11+02 

-♦ 8075+03 

• 325 P +02 

-•0210+03 

*3152+02 

-.0116+03 
• 3309+32 

-.8571+03 

♦3116+02 

-♦8700+03 

*3511+02 

-.0823+33 

•3429+02 

-.8936+03 

.3710+02 

-•9012+03 
. 3016+02 

-.91 36+03 
•3915+02 

-.9225+03 

.<*236+02 

V328 + C3 
•HI IS+Q2 


• V 1 83 + *'2 


MART 

HR 

TE 

BTU/Lft*) • 


<D€4 R 

Qt 

RE LAM 

RETOR* 

r £ci ( 

r2r otHt * s . 

J^LtSS 

P. l**4+03 

Aj*.S»3*0^ 

Jl^P.1061+0! 



gj. 1700+02 

-•$097+03 

-.5553+01 

•3919+01 

•5*66+01 

•1014+02 

*1374+02 

-.52*1*01 

-•5553+01 

•37S6+Q! 

•7373+01 

.7703+02 

• *376 + 02 

-•5147+03 

- .5553+01 

•3417+01 

♦ 9 159+0 K 

•12*5+03 

• 1 l0*+03 

••5626+03 

-.5553+01 

•3199+01 

• 1 1 1 <1*02 

•1961+03 

•2187+03 

-•5766+03 

-.5553+01 

•3346+01 

•l32|+02 

•2937+03 

•3206+03 

• 5*80 + 03 

-.5553+0! 

•3251+0! 

• 1532+0 2 

.1061+03 

•1195+03 

■•5997+03 

-.5553+01 

•3110*0! 

• 1713 + 02 

.5125+03 

•4060+03 

>.6006 + 03 

-.5570+01 

♦ 3054 + 0 1 

•1925*02 

.7013+03 

•7960+0* 

>•6173 + 03 

-.5573+01 

•2973+pl 

• 2 1 3 2 + 0 2 

,0*61+03 

• 1012+01 

■ • 625 1 +03 

-.S57Q+0! 

« 2097+01 

• 2313*02 

.1091+0! 

•1259+01 

.6321+03 

-.5570+0! 

•2026+01 

•2556+02 

. 1 327+0! 

• 1511+Ql 

> « 6 39 1 +g3 

-.5570+C! 

.2761+01 

.2749+02 

• 1587 + 01 

• 1856 + 0! 

•6153+C3 

-.557C+0! 

• 2700+01 

•2971+02 

. 1871 + 0! 

• 22 1 0*0! 

• 6 5 l | + P 3 

-,5570+ri 

.26!! + fl! 

*3183* C 2 

.2188+?! 

•2599+nl 

• 656 q + 03 

5570+01 

•25! |+C1 

• 33d 1 +02 

• 2528+T 1 

•30?3+r! 


LMSC-HREC TM D3 90450 



LMSC-HREC TM D390450 


These are the control variables used in the inviscid surface 
streamline tracing section of the program. These data are 
contained in card 1 of the input. 

These are the initial conditions required to start the tracing 
of each streamline utilizing the fourth order Runge-Kutta 
method. These data are contained in cards 2 and 3 of the 
input. 

These are the control variables used in the convective heating 
rate section of the program. These data are contained in 
cards 1, 4 and 5 of the input. 

The output of the program is identified by the symbols ( 5 ) 
through (25) . These are the body local flow properties, 

heating rates, and characteristic dimensions calculated along 
each inviscid surface streamline. The variables are defined 
a.s follows: 

Indicates the number of the inviscid surface streamline 
presently being traced. 

Axial location on the body surface 

Angular location on the body surface 

Local velocity 

Local entropy 

Local Mach number 

Local density 

Laminar adiabatic wall enthalpy 
Turbulent adiabatic wall enthalpy 
Local wall enthalpy 

Local flowfield temperature at the edge of the boundary layer 
Local impingement pressure 
Local viscosity 
Local Prandtl number 
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Local metric coefficient 
Laminar characteristic running length 
Turbulent characteristic running length 
Laminar convective heating rate 
Turbulent convective heating rate 
Laminar Reynolds number 
Turbulent Reynolds number 

6.3 SAMPLE RESULTS 

Some of the pertinent data contained in Section 6.2 is presented graph- 
ically in this section. The impingement pressure distribution and other local 
flow field properties used by the program were calculated by the Lockheed- 
Huntsville PLIMP program. Figure 6-2 is a plot of inviscid surface stream- 
line coordinates versus X) for several initial streamline starting angles, 

8 . Figure 6-3 is a plot of metric coefficient versus distance X for one of 
these streamlines (0 = .05 radians). In Fig. 6-2, the streamlines diverge at 
an increasing rate as they proceed axially down the cylinder. The rate of 
divergence at any point on the cylinder increases with increasing initial 
starting angle, 0. This divergence, i.e., the spreading of the streamlines, 
is reflected by the metric coefficient plot in Fig. 6-3. Figure 6-4 is a plot 
of the laminar impingement heating rate distribution along the streamline 
with an initial starting angle, 0 = .05 radians. 
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Section 7 

CONCLUSIONS AND RECOMMENDATIONS 

The Streamline Divergence program provides a means to make: (1) 
comparative evaluations of plume impingement due to different propulsion 
system configurations; (2) parametric studies; and (3) vehicle design studies. 
The Streamline Divergence program has demonstrated the capability to trace 
inviscid surface streamlines and calculate "outflow corrected" convective 
heating rates (laminar and turbulent) on surfaces subjected to exhaust plume 
impingement (or any other supersonic flow field). F.R. Dejarnette's stream- 
line divergence methods are applied to determine the streamline patterns and 
calculate the corresponding "equivalent radii" or metric coefficients along 
each streamline for the particular geometry subjected to exhaust plume im- 
pingement. For purposes of demonstration the streamline divergence method 
has been applied to that of exhaust plume impingement on a cylinder. The 
incorporation of other geometries into the computer code is a straightforward 
procedure. 

Convective heating rates are obtained from the boundary layer analysis 
for both the laminar and turbulent flow assumptions along each streamline. 
Eckert's reference enthalpy method is used to account for the temperature 
dependent property variations through the boundary layer. The local Reynolds 
number is output in order for the program user to choose the most applicable 
of the two heating rates. 

Eckert's reference enthalpy menthod has been chosen for use in the 
convective heating rate analysis for demonstration purposes. Other methods 
do exist that may be applied to the problem of exhaust plume impingement 
and, can be incorporated into the computer code with a minimum of program- 
ming difficulty. 
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